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ABSTRACT 


The purpose of this research is to explore the need for 
time-dependent sea-surface temperatures in atmospheric model 
predictions to 10 days. Six and nine-layer versions of the 
Navy Operational Global Prediction System (NOGAPS) are used 
in this study. Control forecasts were made in which the 
_sea-surface temperature (SST) is fixed in time. Test hind- 
casts were made in which the SST was updated at each time 
step of the atmospheric model using interpolations of 
12-hourly SST analyses. The 10-day predictions are compared 
to determine any improvement or degradation due to the time- 
dependent SST. Two cases are analyzed, one during November 
1983 and another duringpwiperil i1964- Use of the tite- 
dependent SST's resuited in significant changes in the fore- 
cast fields of surface heat fluxes and precipitation which 
were phySicaily consistent with the SST trend. Anaiysis of 
15 storm forecasts revealed Significant changes of storm 
track, duration or cyclogenesis in only 4 cases. Three of 
these cases were forecast by the nine-layer version of 
NOSAPS during the April period and one case was forecast by 


the Six-layer NOGAPS during the November period. 
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I. INTRODUCTION 


Atmosphere and ocean interaction 1S recognized as an 
important component of the atmospheric circulation on time 
scales of a month and longer. Air-sea interactions on 
Shorter time scales, such as those over which numerical 
weather prediction (NWP) models are integrated, are less 
well understood. Operational atmospheric prediction models 
generally employ one-way ocean to air influence derived fron 
a time-independent sea-surface temperature (SST) field. The 
quality of forecasts from such models declines rapidly when 
extended beyond five or six days. Many factors contribute 
to this observed degradation: inaccurate specification of 
the initial conditions; inadeguate horizontal and vertical 
resolutions; inaccurate or simplistic parameterizations of 

atmospheric processes; and lack of two-way interaction 
"between the atmosphere and ocean. 

Sea-~surface temperature changes alter the bounds. 
forcing of atmospheric prediction models through the surface 
fluxes of sensible heat, latent heat and momentum. The heat 
fluxes are functions of the air-sea temperature difference 
and planetary boundary layer (PBL) conditions. These fluxes 
can also affect the radiation balance through parameteriza- 
tions of cloud cover and water vapor content. The magnitude 
of the atmospheric response to changing SST's should be 
related to the magnitude and locations of SST changes. The 
nonlinear character of these relationships reguires the use 
of a glokal atmospheric prediction model with sophisticated 
parameterizations of the)? aceracerons between surface 
fiuxes, clouds, precipitation and radiative fluxes. 

This work 1s the second in a series of case studies 


designed to study the necessity and feasibility of coupling 
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atmospheric and oceanic models. The first case study 
(Ranelli, 1984) examined the response of 10-day atmospheric 
forecasts to observed SST changes introduced through sSinmu- 
lated non-synchronous coupling. The temperatures at the sea 
surface were updated by replacing the initial SST's with the 
Fleet Numerical Oceanography Center (FNOC) SST analyses 
valid at 12-hour intervals during the forecast period. In 
this second step, Simulated Synchronous coupling is used to 
introduce observed SST's into the atmospheric modei. Two 
cases, one each in the spring and autumn transition periods, 
are investigated. While a global model 1S used, the anal- 
ysis iS restricted to the northern hemisphere mid-latitude 
oceanic areas and the effects of changing SST‘'s on the 


atmospheric model 10-day predictions. 


3 


Il. BACKGROUND 

Current operational NWP models fall far short o£ the 
approximately 15-day theoretical limit 9f atmospheric 
predictability (Rosmond et al. 19a. Many methods are 
being proposed and tested to incrementally improve existing 
moiels. These include improved numerical differencing 
schemes, spectral model formulation, increasei vertical and 
horizontal resolution, improved initialization techniques, 
and more accurate parameterization of diabatic processes. 
While errforts in these and other areas will ultimately be 
reguired inthe process of NWP model improvement, the 
present study will concentrate on the effects of diabatic 
processes in NWP nodels. 

Diabatic processes become more important in determining 
atmospheric circulations at long time scales. At shorter 
time scales (1-3 days), dynamic processes seem adeguate as 
the primary forcing in NWP models. Diabatic processes begin 
to become important in medium time scales (5-15 days), and 
tend to dominate the long time scale (15-30 days and beyond) 
forecast problen. Pecause the ocean i1isS a major heat and 
moisture source, the air-sea fluxes are very important in 
modeling diabatic processes. 

Latent and sensible heat fluxes across the sea surface 
are modelled as functions of air-sea temperature differ- 
ences. While Sandgathe (1981) concluded that maritime 
cyclogenesis reguired accurate specification of air-sea 
fluxes, no operational model includes time-dependent SST. 
The European Center for Medium-range Weather Forecasts 
(ECMEF) and the United States National Meteorological Center 
(NMC) use the NMC weekiy global SST analysis in initializa- 
tion. The FNOC NOGAPS model is initialized by 12-hourly SST 
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analyses produced by the Thermodynamic Ocean Prediction and 
Expanded Ocean Thermal Structure Systems GLGPS— EO £5) 
described by Clancy and Pollack (1983). Without a time- 
dependent specification of SST, the model representation of 
air-sea interface fluxes and cyclogenesis is bound to dete- 
riorate with increasing forecast time. 

Moiiticar: wal changes yscan™ ocemrm during 5-15 day 
periods. Atmospheric forcing causes large SST changes 
during the spring and autumn transition periods (Camp and 
Elsberry, 1978; Elsberry and Camp, 1978; Elsberry and Raney, 
73). During these periods, warm shallow oceanic mixed 
layers are mechanically mixed, deepen and cool. Winter SST 
changes are generally smaller than transition changes 
because of a deeper ocean mixed layer. Summer SST changes 
are Smaller due to decreased mechanical mixing (fewer and 
less intense storms). The intensity and orientation cf SST 
gradients associated with western boundary currents are 
important in determining cyclogenesis in the western ocean 
basins (Sanders and Gyakun, i Shee Sign riaceant — So. 
increases in tropicalfequatorial regions are manifest in 
local enhancements in the convective cloud amounts on short 
(1-2 days) time scales. 

The specification of a time-varying SST to a NWP model 
would require a corresponding oceanic prediction model. The 
inclusion of atmospheric forcing influences on SST changes 
would reguire an input to the oceanic model from the atmos- 
pheric model. The sensitivity of both models to the input 
from the other must be well understood before truly coupled 


atmospheric-oceanic models become operational. 
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A. COUPLING SCHEMES 


Three methods of model coupling were described by 


Elsberry et al. (19 eZ yr 


1. minimal feedback; 
2. non-SyYRchPoOnous "empl i Wes ane 


3. synchronous coupling. 


Minimal feedback (Fig. la) ais the method used by current 
operational models. It involves specification of SST at the 
initial time and no time-dependence of the SST field. 
Non-synchronous coupling (Fig. 1b) involves SST input fron 
an . essentially independent ocean pDisedictiom model. 
Synchronous coupling (Fig. ic) involves concurrently ruUuieae 
oceahic and atmospheric models which provide boundary 
forcing to each other. 

Ah understanding of the effects of more sophisticated 
coupling schemes on NWP modelS is reyuired before these 
schemes are implemented in operational models. The addi- 
tional degrees of freedom afforded by the new time dependen- 
cieS in bourdary forcing can possibly degrade model forecast 
Skiil, especially that achieved through "tuning" of parame- 
terized processes. Atmospheric model biasSes aay result from 
biases in the oceanic prediction models or improper coupling 
between models. A large number of model verification 
Studies will be required to identify these biases with any 


Statistically significant degree of certainty. 


Be PREVIOUS STUDIES 


Many previous studies have concentrated on long time 
scale (climatalogical) effects of SST changes on the atnos- 
phere. A review of these studies is contained in Elsberry 
et al. (1982). 
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Arpe (1681) examined the effect of artificial SST anonma- 
lies on the performance of the ECMWF model. He showed model 
sensitivity at forecast times of six days and increased 
cyclogenesis with large-scale positive SSI anomalies. 
Ranelli (1984) examined the effects of observed SST devia- 
tions (with respect to the initial conditions) on the NOGAPS 
model. Clear differences in model behavior were noticed 
over a 10-day forecast. The model employing the’ tinme- 
dependent SST seemed to show improved skill in predicting 


mid-latitude maritime storm lifetimes. 


ey 


Iii. PROCEDURE 


The experimental procedure was similar to that of 


Ranelli (1984). Two versions of NOGAPS were used in Ciege 
research. In each of two cases, two NOGAPS integrations to 
10 days forecast time were made. The first of the integra- 


tions was the control run in which the SST was held 
constant. In the second integration, designated the SST 
run, global SST‘'s were updated at every time step (4 
minutes) of the model. The model atmospheric responses of 
each forecast were analyzed for differences. The analysis 


was restricted to the following fields: 


a. surface pressure; 

be. sea-Surifiace temperature; 

ce. surface sensible heat flux; 

d. surface latent heat flux; 

e. surface total heat flux (includes radiation) ; 
£. Convective prece1 cite cmon ana 


ge large-scale precipitation. 


The NOGAPS model was made available by Dr. T. Rosmond of 


the Naval Environmental Prediction Research Facility 
(NEPRF). The general characteristics of the model were 
descrited by Ranelli (1984). During the course of this 


experiment, major changes were made to NOGAPS. The November 
1983 case study waS run with a six-layer version of NOGA?PS 
very Similar to that used by Ranelli (1984). The Apriiiegs 
case was run with the new nine-layer version of NOGAPS. In 
addition to vertical resolution, the two model versions 
differed in the method by which diabatic heating was 
applied. In the earlier version, the diabatic processes 


were calculated and applied once every 10 model time steps 
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(40 minutes). 


The new version also calculates the diabatics 


every 10 model time steps but applies this heating uniformly 


at each time step. 
be applied smoothly with less 
in which diabatic effects are 
ffects of time-dependent SST 


smoothly in the new NOGAPS. 


Conseguently, 


the diabatic effects will 
Shock between the time steps 
The 


introduced more 


ama are Not introduced. 
will also be 


Since additional model levels 


are included, the coupling between the planetary boundary 
layer and the free atmosphere should also be enhanced. 
Cbserved sea-surface temperatures were useji to simulate 
the input from an ocean prediction model. The SST analyses 
of 0000 GMT and 
1200 GMT. to create fields of 
12-hour forward difference SST change. NOGAPS 


difference SST 


are produced twice daily for valid times 


These analyses were used 
was modified 
to input these 12-hour forward Changes and 
SST change at each time 


Ou 


apply an incremental portion of the 


step. This process simulates a synchronous coupling 


NOGAPS to a "perfect prog" ocean model. 


The NOGAPS initial conditions used for tne case studies 


were identical to those used for the operational forecasts 


with the same starting times. The case starting times were: 


ae 1200 GMT 31 October 1983; and 
De 1200) GNI 21 April 1984. 
The model output parameters were analyzed by three 
methods. Forecast and analyzed vortices were objectively 


tracked using the Systematic Errors Identification System 
(SEIS) developed by Brody, et al. (loco) at NEPRPS SEIS is 
basead on the vortex representation system of Williamson 
(1981). The system tracks vortices in a 900 grid point 
‘window’ ain the FNOC Co oom POOlareestcreOgrLadphic grid. 
Separate SI runs were reguired to track Pacific and 
Atlantic storms. The other model output fields were 
analyzed in Atlantic and Pacific ‘windows' of the FNOC 


eo 


73x144 spherical grid which extended from 20N to 60N. The 
SST and control run forecast parameters were differenced 
(SST-control) and both contour plots and field statistics 
were used to determine the relationships between the fore- 


cast differences and the time-dependent SST. 
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IV. HEAT FLUX ANALYSIS 

The introduction of time-dependent sea-surface tempera- 
tures into the NOGAPS model results in immediate changes in 
the surface fluxes of sensible and latent heat. These 
fluxes in NOGAPS are parameterized according to Deardorff 
mo 7 2). Because the PBL and convective parameterization 
schemes in NOGAPS are intimately coupled, the changes in 
air-sea fluxes may change the diabatic forcing throughout 


the troposphere. 


Ae SEA-SURFACE TEMPERATURE CHANGES 


In both the November and April cases, Signe cant SST 
changes evolved during the 10-day forecast period. mue SST 
difference field (day 10- initial conditions) for the 
Wovember case iS Shown in Fig. 2. Large-scale negative SST 
changes are evident in the western Pacific Ocean, Sea of 
Okhotsk, eastern Pacific Ocean, and mid-Atlantic Ocean. 
There were small areas of weak positive SST change to the 
west of the British Isles and to the east of the north- 
eastern United States. $$The SST difference fields for the 
April case are Shown in Fig. Se While both positive and 
negative SST changes are shown, the overall mean SST differ- 
ence if poSitive. Large positive differences are found in 
the Sea of Japan, Yellow Sea, eastern Pacific and western 
Atlantic Oceans. Negative differences (decreasing SST with 
time) occur in the western Pacific (east of Japan) and’ north 
Atlantic (south of Greenland). The evolution of the field 
mean SST differences for each ocean basin in both cases is 
Shown in Fig. 4. The 10-day trends in the mean SST differ- 


ences are essentially linear. Correlation coefficients of 
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linear regressions on the field-mean SST differences were 
between 0.86 and 0.99. Some of this trend is due to the 
actual SST observations. JE aljelfalah jest (ey8) 7 the TOPS-EGa#s 
forecast/analysis scheme includes a time-dependent effect in 
the first guess and a reverSion to climatology in data- 
Sparse areas. Each of these effects will also contribute to 


a general trend. 


Be. SURFACE HEAT FLUX CHANGES 


As in Fanelli (1984), the surface heat fluxes were 
responsive to the SST changes. Howeve.c, the spatial 
patterns of heat flux differences were not always related 
directly to the patterns of SST difference. The heat flux 
difference patterns often took the form of couplets of posi- 


tive and negative differences. 
i. Sensible Heat Flux 


The field mean heat flux differences did show trends 
SiMiiar to the SST differences. The 12-hour fields of 
sensible heat flux differences were extremely noisy. 
Features in the difference fields were rarely stable for 48 
hours. The largest differences, both positive and negative, 
were observed in the western ocean basins in the November 
case and in the Pacific Ocean south of 30 N for the April 
case. 

The development of the mean sensible heat flux 
differences iS shown in Figs. 5 through 8. While the field- 
mean diitferences reach 5-20 percent of the control forecast 
mean sensible heat flux, peak values in the difference 
fields exceed 50 percent of peak values in the control fore- 
cast fields. Peak differences exceed the mean control 
fluxes by factors “of Zao. Correlation coefficients of 
linear regressions of the field-mean differences were irom 
Os Otte. 0.567 


2. Latent Heat Fiux 


The latent heat flux difference fields showed less 
noise and more stability than those of sensible heat flux. 
A direct correlation of the latent flux difference fields 
and SST difference fields was not obvious. Representative 
latent heat flux difference fields are Shown in Fig. 9. In 
the November case, strong poSitive-negative couplets devel- 
oped in the western Atlantic region. The Pacific Ocean 
areas south of 30 N exhibited the largest differences in the 
Aprii case. Analysis of the field statistics (Figs. 10 
through 13) shows mean differences reaching 5-12 percent of 
the control means. Peak values inthe difference fields 
reach 25-75 percent of the peak values of the corresponding 
control latent heat flux fields and exceed fieli-mean fluxes 
or timemeonteol run by £actors of 2-6. Linear regression of 
the field-means produced correlation coefficients ot 
0.75-0.88 in the November case and only 0.22-0.60 in the 
April case. The differences in latent heat fluxes after 10 
days during November are at least 10 times as large as the 
dizferences in sensible heat fluxes, and represent consider- 


able reductions in evaporation over the entire field. 
3. Total Heat Flux 


NOGAPS total heat flux includes the sensible and 
datent fluxes plus the net radiative flux at the surface. 
The radiation is affected by model cloudiness. Fields of 
total heat flux difference tend to be correlated with fields 
of latent heat flux aqifference. Field statistics (Figs. 14 
through 17) showed mean differences of 8-20 percent over the 
10-day forecast periods. Peak differences at specific tines 
were 40-100 percent of the control peak fluxes and exceeded 
control field-means by factors of 2-10. Linear regression 


of the field-means produced correlation coefficients of 
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O0.44-0.94 in the November case and 0.12-0.46 in the April 
case. Especially in the November case, these differences 
after 10 days represent Significant changes in the surface 
heat budget terms which wouid be used to drive the ocean 


model. 


C. SUAMARY 


The introduction of time-dependent sea-surface tempera- 
tures resulted in significant field-mean sensible, latent 
and total heat flux differences of up to 20 percent over 
10-day forecast periods. Maxima in the flux differences 
between the SST and control runs after 10 days were of the 


Same order of magnitude as peak values in the control fiux 


fields. Thus, the surface fluxes at individual points may 
differ significantly from _the scontrol run when time- 
dependent SSI's are specified. Stable features in the 


difference fields of latent and total heat flux were 
observed to intensify and move consistently during the 
second half of the forecast period. The largest changes 
were noted in the western Atlantic in the November case and 
in the Pacific (south of 30 N) for the April case. 

The forecast differences in heat fluxes are of such 
magnitude that they can be expected to change the nature of 
the atmospheric (and oceanic) forecast at some future time. 
Local changes in the heat fluxes could also be expected to 
seriously affect modei-derived forecasts, such as those 
produced by model-output statistics and systems Similar to 
the Navy Operational Local Atmospheric Prediction System 
(NOLAPS), Since these products depend on the predicted 
surrace fluxes. | 

The changes in the surface heat fluxes were greater than 
those observed by Ranelli (1984). Three key differences in 


the experiment procedures account for the differences in 
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heat flux observations. 


Peete 6 6CayS OL SST 


input vice 7 days during the Ranelli study. The 


of diabatic heating by the new NOSAPS increased 


or a time-dependent 


minutes 


(vice 12 hours) 


Sodas didethe updating of 


of forecast time. 
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Change were 
application 
the effects 
SST every 4 


V. PRECIPITATION ANALYSIS 

Changes in the surface heat fluxes induced by time- 
dependent sea-surface temperatures can be expected to affect 
convective and large-scale precipitation forecasts. A study 
by Bosse (1984) Showed that surface neat fluxes could be 
related to NOGAPS large-scaie precipitation fields through 


layer cioud instability. 


A. CUBULUS PRECIPITATION 


The development of the cumulus precipitation difference 
fields during the 10-day forecasts 1S shown int Figs. 12 
through 21. Over both the North Pacacic and North At lami 
there are 25'to 50 percent decreases in mean cumulus precip- 
itation during the November forecast. In the April fore- 
cast, mean cumulus precipitation over the Pacific decreased 
almost 15 percent while cumuluS precipitation over the 
Atlantic increased 20 percent. Peak diitference vaiues 
exceeded 100 percent of peak at some time during each of the 
forecasts. The decrease in cumulus precipitation during the 
November period seems logical in view of the decreasing mean 
sea-surface temperatures and latent heat fluxes. However, 
the magnitude of the “changes as SUcprisiagdy Danger 
Correlation coefficients of 0.74 were obtained in linear 
regression of both the Atlantic and Pacific mean cumulus 
precipitation differences. The magnitudes of the SST and 
latent heat flux field-mean differences were much less in 
the April case than in the November case. The variance of 
both mean SST change and mean flux change was greater in the 
April case than in the November case due to the presence of 


areaS of strong positive and negative SST change in both 
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ocean basins. This may explain the less consistent changes 
in cumulus PECCIPitabioen saquring the April forecast. 
Correlation coefficients for the linear regression of the 
mean cumulus differences during the April period were 0.61 


meeeti@e Pacitiac and 0-007 in the Atlantic. 


Be. LAKGE-SCALE PRECIPITATION 


The aevelopment of the large-scale precipitation differ- 
ence fields iS Shown in Figs. 22 through 25. Mean differ- 
ences of large-scale precipitation are less than 10 percent 
except ror the Aprii Atlantic at the end of the forecast 
period. Peak differences range from 75 to 100 percent of 


peak control values. 


G  UBMARY 


Significant mean and peak precipitation differences were 
observed during the 10-day forecasts. AS with the heat flux 
differences, these Changes in modei precipitation grow with 
time. Long range model-derived products such as nodel- 


Output statistics wiil be affected by these differences. 
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VI. ANALYSIS OF STORM TRACKS 

The SEIS program (Brody, et al., 1984) was uSed to oODjieg= 
tively analyze the position, size, Shape and amplitude of 
low pressure systems in the control and SST ‘forecasts and 
the analyses. The analysis was performed on surface pres- 
sure only. 

the SEIS storm parameters are defined in a tanner 
Similar to that of Willianson era). Amplitude is the 
difference between the storm central pressure and the 
latitudinally-averaged surface pressure. Radius is the size 
of the storm expressed in terms of the spacing of the FNOC 
Northern Hemisphere polar stereographic grid. ELLIpt esate 
1s the square root of the ratio of the Semi-major and semi- 
minor axes of the ellipse which best fits the analyzed storn 
Shape. Angle is the angle between th2 x axis and the sStorn 
ellipse major axis, measured counterclockwise from the x 
axis. 

in the SEIS program, the forecast low pressure systems 
are analyzed first, then low pressure systems in the anal- 
yses, then forecast and analyzed systems are matched 
according to set criteria. The nominal SEIS matching 
criteria are plus or minus 2.5 grid points in both the x and 
y directions and plus or minus 12 hours in time. SERSwars 
routinely used to evaluate 60 hour forecasts and the 
matching criteria are optimized for this period. AS a 
result, the vortex tracking and description portions of SEIS 
perform satisfactorily throughout the 10-day forecastS, (jaa 
the automatic matching portion of SEIS becomes ineflecuame 
after five days. This is primarily a function of the ove 
restrictive time matching criteria. inspection of the SEIS 


output indicated that several gooi storm track matches in 
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the 5-10 day portion of the forecast period could be made if 
the time matching criteria were increased to plus or Minus 
48 hours. Manual matching with the standard distance 
criteria and 48 hour time criteria was used in the 5-10 day 
portion of the forecasts. 

Labels were aSsSigned to’ the matched storm tracks for 
discussion and reference purposes. The storm labels consist 
eo Chdlacters. The first Character is the first letter of 
the forecast month (N ror November, A for April); the second 
character is the first letter of the ocean basin (A_ for 
merantic, © Lor Pacific); and the third character isa 
humeral assigned on the basis of chronological order of 
storm identification within a particular forecast and ocean 


baSin. 


A. STORMS DURING NOVEMBER 


The November control and SST forecasts were made with 
the Ssix-level version of NOGAPS. A total of 15 storms with 
lifetimes greater than 36 hours were tracked in the surface 
pressure analyses. Fifteen stormS were present in the 
control forecast while 16 were tracked in tne SST forecast. 
Matches were obtained for eight of the 15 analyzed storm 
packs. 

Relatively minor differences between the control and SST 
torecast storm positions and parameters were noted. These 
qifrferences grew with time. The tracks for storm NA5, which 
showed the largest track and parameter differences, are 
Shown in Fig. 26. This’ storm was first identified at 72 
hours in the analyses~ and at 96 hours in both forecasts. 
The analysis and control forecast storm lifetimes of 120 
hours are identical whiie the SST storm lifetime was 24 
hours longer. The forecast storm tracks diverge after the 


third forecast but remain almost parallel. The control run 
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storm remains quaSi-stationary north of Great Britain four 
days aiter first detection, while the SST run storm stalls 
for only 12 hours in the same region and taen proceeds to 
the northeast. The cyclogenesis and early storm tracks of 
NA5 are in the region of maximum SST change, but the late 


tracks and cyclolysis positions are in the region of mininun 


SSimchaniges The SEIS parameters for storm NAS are shown in 
jae (aera The forecast amplitude and period of maximum 
radius lag the analyzed parameters by 48 to 60 hours. The 


control run waS more accurate than the SST run in this 
particular case, although the aifferences are much smaller 
than the departures of the forecasts from the analyses. The 
NA5 forecast was the only one of the eight November storms 
for which one forecast could be judged more accurate than 
tie other. Tae descriptions, tracks and parameters for the 


other seven November storms are contained in Appendix A. 


Be. STORMS DURING APRIL 


The April control and SST forecasts were made with the 
hnine-layer version of NOGAPS. A total of 11 storms with 
lifetimes greater than or egual to 36 hours were tracked in 
the surface pressure analyses. Fourteen were tracked in the 
control forecasts, while 15 were present in the SST fore- 
casts. Matches were obtained for seven analyzed storm 
tracks. Significant differences between control and SST 
forecast tracks and parameters were found for three of the 
seven matched tracks. Those matched tracks with significant 
differences are discussed here. The remainder of the April 


storm tracks are discussed in Appendix A. 


1. Storm AA3 


SS Se Se SS ESS 


The tracks and SEIS parameters for storm AA3 are 


Showh iD Fags: 26 qaniemeee This storm was first detected in 
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the analyses at 84 hours and was forecast only by the SST 
Fume at’ 108 hours. The forecast Storm was weaker (5 mb) and 
shorter-lived (36 vice 168 hours) than the actual storm, but 
the cyclogenesis position and initial track were well 
predicted. This was the only case for which only one of the 
forecast models predicted a matched storm track when the 
other forecast model did not. The cyclogenesis position was 
to the west of, and the storm track over, an area of maxinun 
positive SST change. The SST forecast was clearly superior 
to the control run in this case, especially considering that 


the storm developed around the fourth day of the forecast. 


2. Storm AA4Y 


The tracks and SEIS parameters for storm AA4Y are 
SiowWh in Figs. 30 and 31. Storm AAY was forecast by both 
the control and SST runs at 120 hours and was first detected 
in the analysis 12 hours later. The cyclogenesis position 
waS in an area of maximum positive SST change (1.2 C). Both 
forecasts overestimated the actual storm Jlifetime of 72 
hours by 48 hours. The forecast tracks were Similar until 
the sixth forecast (third day after formation), when the SST 
storm began erratic movement. The forecast SEIS parameters 
ier 19. 31 are Similar, with significant differences in 
radius and angle at the second forecast period and radius 
and ellipticity at the eighth forecast period. The control 
track forecast for AA4 was very slightly more accurate than 
Grew ou track Lompecast during the first three days. The 
differences between the forecastS again are much less than 
the departures of each forecast from the analysis. 


3. Storm AP3 


SS ee ce 


The tracks and SEIS parameters for storm AP3 are 
shown in Figs. 32 and 33. Cyclogenesis occurred at 60 hours 


in both forecasts andthe analysis. The cvclogenesis and 
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tracks were 1n alakea of poSatave SST Changess(Om4-0- cae. 
The control run forecast a weak Storm with a lifetime of 
orly 12 hours whereas the SST run predicted a lifetime of 60 
hours. Since the analyzed storm lifetime was 48 hours, the 
SST run more accurately predicted the storm lifetime and 
track in this case. The amplitude (Fig. ~33) of thisenediee 
tively weak storm was forecast rather accurately in the SST 


rune 


C. SUMMARY 


Significant differences in the forecast tracks and SEIS 
parameters were noted in four of the 15 matched storn 
tracks. The control run waS more accurate in two of the 
cases, correctly forecasting November cyclolysis when the 
SST run did not and more closely matching an April storm 
track when the SST run exhibited erratic movement. The SST 
run waS more accurate in two April cases, forecasting an 
analyzed cyclogenesis which the control run did not and 
Maintaining a storm in which the control run forecast prema- 
ture cyclolysis. The different method of applying the 
diabatic forcing andthe increased vertical resolution of 
the nine-level NOGAPS used in the April runs made direct 
comparison of the November and April forecasts difficult. 
The April resultS were very encouraging. it appeared that 
the nine-level NOGAPS was sensitive to the effects of the 
time-dependent SST and that this sensitivity translated into 


more accurate cyclone predictions in two cases. 
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VII. CONCLUSIONS 
The introduction of time-dependent sea-surface tempera- 
tures into the NOGAPS model resulted in significant changes 
in mean and peak surface heat fluxes and precipitation. 
Forecasts which use these fields directly, such as model- 
output statistics, TOPS and NOLAPS-type models, may benefit 
fron the incorporation of time-dependent sea-surface 


mem peratures. 


Differences between the control and SST forecast storn. 


tracks and SEIS parameters were significant in 4 of the 15 
matched storms. The April SST run, which utiiized the nine- 
level NOGAPS model, correctly predicted cyclogenesis' and 
cyclone maintenance for two cases in which the control run 
ia not. The control run more accurately predicted a 
November cyclolysis event and one April storm track. Only 
minor differences in forecast positions and SEIS paraneters 
were noted in the other 11 matched storms. Differences in 
the SEIS parameters were evident before those in storm posi- 
fron . The new nine-level NOGAPS, with sophisticatic appli- 
cation of diabatic forcing, seems sensitive to the effects 
of time-dependent SST in away which improves the accuracy 
of some forecasts. 

This study was limited in its application and results. 
Further studies on the response of atmospheric prediction 
woielis to time-dependent SST will be reguired to make 
statistically Sigmiticant inferences about the effects of 
such boundary conditions on the guality of the numerical 
forecasts. Additional case studies with the analyses 
extended to cover the tropical and Southern Hemisphere ocean 
areas should be conducted. Additional model variables, such 


as cloudiness, PBL thickness and static energy and diabatic 
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and cumulus heating should be analyzed. Further attempts 
Should be made to verify the heat flux and precipitation 
forecasts of models which incorporate time-dependent SST to 
determine model sensitivity and biases. The effect of model 
output changes on proposed long range model-output statis- 
tics and other model-derived forecasts should also be 


evaluated. 
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APPENDIX A 
STORM TRACK DATA 


The 11 forecast storm tracks which were judged to show 
jess Significant differences are discussed in this appendix. 
in general, storms which were present in the forecast 
initial conditions or developed early in the forecast period 
did not develop significant differences between the control 


and SST forecasts. 


A. STORHS DURING NOVEMBER 


1. Storm NAT 


The tracks and SEIS parameters for storm NA1 are 
shown in Figs. 34 anG 35. This storm was present in the 
initial conditions of both forecasts. Only  Oinon d2ztter— 
ences between the SST and control forecast parameters 
develop after six to eight 12-hour forecasts. Both forecast 
storms moved too slowly after 24 hours and persisted longer 


than the analyzed storn. 


2. ~eOumeslie2 


The tracks and parameters for storm NA2 are Shown in 
Pegs. 36 vand 37. In this case both forecasts accurately 
predicted storm lifetime while moving the storm too quickly 
and in a more zonal path than was analyzed. Menon ditt er- 
ences in only two parameters, position and angle, can be 


seen at the sixth forecast period. 


oie Storm NA3 


The tracks and parameters for storm NA3 are shown in 


pigs. soeand, 39. This storm formed at 72 hours into the 
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Z240-hour forecasts. Minor differences in the forecast storm 
parameters are evident from the time of cyclogenesis. Both 
forecasts underestimated storm lifetime. The analyzed storm 
moved erratically and regenerated after the tenth forecast 
period. Neither model was able to predict the  storn 


regeneration. 


4. Storm NAY 


The tracks and parameters for storm NAY are Shown in 
Figs. 40 and 41. This storm formed at 72 hours. While the 
control and SST forecast parameters are Similar, the control 
Storm lifetime was 36 hours greater than that of the SST 
torecast. The analysis and both forecast tracks were 


CEE at 1c- 


De. SEODMmWe a 


The tracks and parameters for storm NP1 are shown in 
Figs. 42 and 43. This Aleutian low was present in the 
initial condition of both forecasts and persisted for 7 days 
in each case. A rather extreme southeastward deflection in 
the analysis occurred between the third and fourth from the 
end positions. Differences in the SEIS forecast parameters 
are apparent after the seventh through ninth forecasts. 


Only minor forecast storm position differences are evident. 


CO.  StOmimenes 


The tracks and parameters for storm NP2 are shown in 


Figs. 44 and 45. This storm formed at 24 hours into the 
240-hour forecast period. Both forecasts underestimate 
storm lifetime. No differences in the forecast storm 


tracks, and only minor differences in the forecast SEIS 


parameters, were noted. 


S0 
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The tracks and parameters of storm NP3 are shown in 
Pato eeooeand 472) Storm NP3 waSethe extratropical transition 
of a tropical cyclone which was not analyzed until 48 hours 
imo the forecast.™ SEIS did not track this vortex until 132 
hours because the storm size was less than 2 grid lengths. 
Both forecasts predict cyclogenesis in the area of the 
g@emlad sextratropical transition anproximately 48 hours 


before it actually occurred. 


Be. STORMS DURING APRIL 


Alize Storm AAT 
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The tracks and parameters for storm AA1l are shown in 
Figs. 48 and 49g. This high Jlatitude storm was present in 
the initial conditions of both forecasts. Despite the 7120 
to 132 hour forecast storm lifetime, only minor differences 
in the forecast positions and parameters are noted. The 
analyzed and forecast storm tracks are roughly parallel to 


the climatic ice edge. 


z.- Storm AA2 


The tracks and parameters for storm AA2 are shown in 
Prose 50 and 51. This storm formed 12 hours into the fore- 
cast period. Both forecasts seriously overestimated the 


analyzed storm ilifetime by of 36 hours by 72 to 84 hours. 


3. Storm API 


The tracks and parameters for storm AP1 are Shown in 
fess oc and 53. This storm formed 24 hours into the fore- 
cast period. Both forecasts moved the storm too far inland 
and overestimated the actual lifetime of 24 hours by 48 
hours. One would not expect significant effects due to 


time-dependent SST on the time scale of these forecasts. 


Si 


G. Storm AP2 


The tracks and parameters for storm AP2 are shown in 
Figs. 54 and 55. This storm formed at 48 hours in the anal- 
ysis and both forecasts. Both forecasts underestimated the 
actual storm lifetime of 168 hours by 96 hours. AS Shown by 
the radius in Fig. 55, this was a rather snall storn. 
Although the radiuS was well predicted for 6-7 12-hour 
periods, the amplitude was too small and the SEIS could no 


longer track the storm. 
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12. Similar to Fig. 10, except during April. 
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14, Field statistics for total heat flux during November in 
the Atlantic region. (a) Field-mean difference in 
m-cal/cn@~h. b) Ratio of field-mean difference to 
ield-mean control. (c) Ratio of peak difference to peak 
control. (d) Ratio of peak difference to mean con t£ous 
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16. Similar to Fig. 14, except for April. 
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19. Similar to Fig. 18, except in the Pacific region. 
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21. Similar to Fig. 20, except in the Pacific region. 
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28. Similar to Fig. 26, except for storm AA3 which 
formed at tau= 84 hours. 
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26, except for storm AA4 which 
formed at tau= 120 
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32. Similar to Fig. 26, except for storm AP3 which 
formed at tau= 60 hours. 
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34. Similar to Fig. 26, except EOResvonmaiat which 
was present in the initial conditions. 
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36. Similar to Fig. 26, except for storm NA2 which 
waS present in the initial conditions. 


o 


74 


AMP MB 


eee) Seas 


Bs 





100 les 200 


90 








Poon eon 
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38. Similar to Fig. 26, except for storm NA3 which 
formed at tau= 72 hours. 
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40. Similar to Fig. 26, except for storm NA4 which 
formed at tau= 72 hours. 
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41. Similar to Fig. 27, except for storm NA. 
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42, Similar to Fig. 26, except for storm NP1 which 
waS present in the initial conditions. 
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uu. Similar to Fig. 26, except for storm NP2 which 
forned at tau= 24 hours. 


82 


0 


ME 


rer 


45. 








Foul 


Stmoear) to Fig. 


150 


1S 20 


lig a 
10 





een 


27, except for storm NP2. 


8 3 





990°o5° nwo° N8S° N40° N4S°N S50 N 


130° "140° £ 150° B 16M ©) W7oeme ume te a/teeeianeen 


46. Similar to Fig. 26, except for storm NP3 which 
formed at tau= 120 hours. 
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Simuganr to Fig. 27, except for storm NP3. 
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Similar to Fig. 26, except for storm AA2 which 
formed at tau= 12 hours. 
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52. Similar to Fig. 26, except for storm AP1 which 
formed at tau= 24 hours. 
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Similar to Fig. 26, except for storm AP2 which 
formed at tau= 48 hours. 
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